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Abstract 

CN . We calculate the partonic cross sections for the hadroproduction of polarized XcJ 

mesons using the factorization formalism of nonrelativistic QCD. We also describe 
■ how the polarization is passed on to the J/ip mesons in the radiative feed-down 

CN . XcJ — > J/ip + r y and in the double-cascade decays of ip' mesons via XcJ intermediate 

\ states. These represent the missing ingredients needed to predict the polarization 

' of prompt J/tp mesons at the Fermilab Tevatron. 

§ ! PACS numbers: 13.85.-t, 13.85.Ni, 14.40.Gx 
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1 Introduction 



Since its discovery in 1974, the J /if) meson has provided a useful laboratory for quanti- 
tative tests of quantum chromodynamics (QCD) and, in particular, of the interplay of 
perturbative and nonperturbative phenomena. The factorization formalism of nonrel- 
ativistic QCD (NRQCD) [jlj provides a rigorous theoretical framework for the descrip- 
tion of heavy-quarkonium production and decay. This formalism implies a separation of 
short-distance coefficients, which can be calculated perturbatively as expansions in the 
strong-coupling constant a s , from long-distance matrix elements (ME's), which must be 
extracted from experiment. The relative importance of the latter can be estimated by 
means of velocity scaling rules, i.e. the ME's are predicted to scale with a definite power of 
the heavy-quark (Q) velocity v in the limit v — > 0. In this way, the theoretical predictions 
are organized as double expansions in a s and v. A crucial feature of this formalism is 
that it takes into account the complete structure of the QQ Fock space, which is spanned 
by the states 2S + 1 L^ with definite spin S, orbital angular momentum L, total angular 
momentum J, and colour multiplicity a — 1, 8. In particular, this formalism predicts the 
existence of colour-octet processes in nature. This means that QQ pairs are produced 
at short distances in colour-octet states and subsequently evolve into physical (colour- 
singlet) quarkonia by the nonperturbative emission of soft gluons. The greatest triumph 
of this formalism was that it was able to correctly describe the cross section of inclu- 
sive charmonium hadropro duct ion measured in pp collisions at the Fermilab Tevatron , 
which had turned out to be more than one order of magnitude in excess of the theoretical 
prediction based on the traditional colour-singlet model (CSM) ||f|]. 

The experimentum crucis of the NRQCD factorization formalism is the measurement 
of the polarization of charmonium with quantum numbers J PC = 1 , i.e. J/if> and 
if)' mesons, in direct hadroproduction, which is predicted to be fully transverse at suffi- 
ciently large values of the transverse momentum (pt) [[J. A convenient measure of the 
polarization is provided by the variable 

a T -2a L 

a = , 1) 



where <tt and <jl are the transverse and longitudinal components of the total cross section 
a = a T + cr L , respectively. The unpolarized case is characterized by a = 0, while a = 1 
and —1 correspond to fully transverse and longitudinal polarizations, respectively. As 
for the J I if> and if)' mesons, a is best measured by analyzing the angular distributions of 
their decays to fi + n~ pairs, which are proportional to 1 + a cos 2 9*, where 9* is the angle 
between the /i + direction in the charmonium rest frame and the charmonium direction in 
the laboratory frame. The first such measurement was recently performed for if)' mesons by 
the CDF Collaboration at the Tevatron 0. Unfortunately, the statistics are still too low 
to allow for a meaningful test of the NRQCD prediction |7|||. In the same experiment, the 
polarization was also measured for prompt J /if) mesons, i.e. those which do not stem from 
the decays of bottom-flavored hadrons, with about a hundred times more statistics. This 
makes it necessary to include in the theoretical prediction the production of polarized J/ if> 



2 



mesons via the decays of the heavier charmonium states Xci, Xc2, and if)', which are directly 
produced. In fact, these feed-down channels make up approximately 15%, 15%, and 10% 
of the prompt- J /if) signal, respectively. There is an additional, yet smaller contribution to 
prompt J /if> production from the double- cascade decays of directly produced if>' mesons 
via Xcj intermediate states. The purpose of this paper is to provide analytic results for 
the missing ingredients which are necessary to make quantitative predictions for all these 
contributions. The resulting prediction for a has already been presented in Ref. ||. It 
is consistent with the CDF data at intermediate values of p?, but it disagrees by about 
three standard deviations from the data point of largest pt- 

In order to predict a for prompt J /if) mesons to lowest order (LO) in the NRQCD 
factorization formalism, we need the polarized cross sections for the partonic subprocesses 
a + b — > cc(n) + d, where a,b,d = g,q,q with q = u,d,s, for n = 3 S[ 1,& \ 3 Pj 1 ' 8 ' ) , and 
1 S , o 8 ' 1 . The corresponding results for the unpolarized case are well established P,|^JT0|. 
The polarized cross sections may be found in Refs. |8|,[lT|, with the exception of those 
for the P-wave colour- singlet cc Fock states n = 3 P±^ and 3 P^ , which are relevant for 
the Xci and Xc2 mesons. These missing results could have already been obtained in the 
CSM, before the advent of the NRQCD factorization formalism. However, the calculation 
is somewhat tedious, especially for the spin-two case of Xc2- The calculation is further 
complicated by the fact that a is measured in the hadronic center-of-mass (CM) frame, 
so that one cannot benefit from the usual simplifications which occur when one chooses 
the partonic CM frame. Instead, the cross section does not only depend on the partonic 
Mandelstam variables, but also on the boost factor relating the partonic and hadronic 
CM frames. Therefore, the helicity amplitude method, which provides a powerful tool 
for the treatment of complicated partonic subprocesses, such as g + g — > Xc2 + g without 
polarization, cannot be directly employed. In Refs. |7],§,|Il[], a Lorentz-covariant formalism 
for describing the polarization of massive spin-one particles in the hadronic CM frame was 
elaborated. We adopt this formalism and extend it to the spin-two case. Furthermore, we 
need the probabilities for the transitions from the cc Fock states of the Xch Xc2, and if>' 
mesons to those of the J /if) meson in the presence of polarization. Owing to heavy-quark 
spin symmetry, they are trivial for the single-cascade decay of the if)' meson. The results 
for the 3 P^ and 3 P^ channels may be found in Ref. |TJ|], where gluon fragmentation 



to polarized Xci and Xc2 mesons with subsequent radiative feed-down to polarized J /if) 
mesons was studied. The missing pieces of information are provided here. 

This paper is organized as follows. In Sec. |2|, we present a general Lorentz-covariant 
formalism which allows us to conveniently describe the partonic cross sections for the 
inclusive production of massive spin-one and spin-two bosons, the helicities of which are 
fixed in the hadronic CM frame. In Sec. |3|, we decompose the cross section for the prompt 
production of longitudinal J /if) mesons into the direct, single- cascade, and double cascade 
contributions, and we explain how the various components factorize into cross sections 
of partonic subprocesses a + b —* cc{n) + d with definite cc helicity, scalar ME's, and 
measured branching fractions. In Sec. |], we present separate theoretical predictions of 
a for J /if> mesons from the single- cascade decays of Xci, Xc2, and if)' mesons, and from 
the double-cascade decays of if)' mesons via Xcj intermediate states. Our conclusions are 
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summarized in Sec. [5| In Appendix we present the ingredients which enter the analytic 
expressions for the polarized cross sections of the partonic subprocesses a + b — > cc{n) + d, 
with n = 3 P^ and 3 P 2 , which represent the missing links for the theoretical prediction 
of a from prompt J /if} mesons. 



2 Covariant Tensor Decomposition of Polarized Cross 
Sections 

The partonic cross sections for the production of massive spin-one and spin-two bosons, 
the polarizations of which are to be measured in the hadronic CM frame, depend not only 
on the partonic Mandelstam variables, but also on the boost factor relating the partonic 
and hadronic CM frames. A Lorentz-covariant formalism to deal with this problem in the 
spin-one case was elaborated in Refs. [J7]||n] . It can be applied to the direct production of 



polarized J/ip, Xa, and if)' mesons. In order to also include Xc2 mesons in the polarization 
analysis, we need to generalize this formalism to the spin- two case. Prior to doing that, 
we recall the spin-one formalism. 

In the following, we work in the hadronic CM frame. Let e M (A) be the polarization 
four-vector of a spin-one boson H with mass M, four-momentum P, and helicity A. In the 
case of longitudinal polarization (A = 0), a Lorentz-covariant expression reads [[7|,p|,[ill : 

'{P ■ Qf - M 2 s 



where Q is the total four- momentum of the colliding hadrons and s = Q 2 . It is convenient 
to decompose the polarization sum, 

1 pfJ. pv 

r = £ e"W(A) = -sr + -^r, (3) 

into longitudinal and transverse components as 

r = E pK> ( 4 ) 

|A|=0 



where 



Specifically, we have 



PfA|= E ^W(A). (5) 

A=±|A| 



fC = z»z v , pT^-pT (6) 

Let M. be the transition matrix element of the partonic subprocess 

a( Pa )+b(p b )^H(P,\) + d, (7) 
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where the four-momenta and helicities given in parentheses refer to the hadronic CM 
frame. If we average the absolute square of M. over the spin and colour states of the 
initial-state partons a and b, and sum over A = ±|A| and the spin and colour states of the 
final-state parton d, we can write the outcome as 



\M 



S, 



where is a rank-two Lorentz tensor which depends on p a , pb, and P. Exploiting the 
the symmetry and orthogonality properties pj^j = pYfi and P^p^ — 0, respectively, we can 
decompose as 



C i S i 7 



(9) 



i=l 



where the common factor F and the coefficients c, are scalar functions of the partonic 
Mandelstam variables 



S = (Pa+Pb)' 
t=( Pa -Pf 

u=(p b - P) 2 



(10) 



which satisfy the relation s + t + u = M 2 , and 



/1U 



6 4 



pZpI 
Pbpb7 
pIpV 



fill 



We now turn to the spin- two case. Let e M ^(A) be the polarization tensor of a spin-two 
boson H with mass M, four-momentum P, and helicity A. It can be constructed from the 
polarization four- vectors of two spin-one bosons with the same mass and four-momentum 



with the help of the addition theorem for two angular momenta as JT3 

i 



(A)= (l,A 1 ;l,A 2 |2,A)e^(A 1 )e^(A 2 ), 

Ai,A 2 =-l 



where (1, Ai; 1, A21 2, A) are Clebsch-Gordon coefficients. Similarly to Eq. 
pose the polarization sum, 



(12) 

), we decom- 



P 



pupa 



e^*(A)e" CT (A) = - (p^p ua + p^p vp ) 



1 



(13) 



as 



P 



pupa 



E pvpa 
P\\\ ' 

|A|=0 



(14) 
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where 



Specifically, we have [12 



f%T = E ^*(A)e^(A). 
A=±|A| 



(15) 



pupa 

Po 



p pupa 



avpa 
H2 



~ (2/C 



pD, 



(p? p pT 



pTp? 



pTpT) 



(16) 



We now consider the case when H in Eq. (^|) represents a spin-two boson. Then, 
Eq. (|SD is replaced by 

(17) 



\M\ 2 



s, 



pupa 
pupaP\\\ , 



where S pvpu is a rank-four Lorentz tensor which depends on p a , pt,, and P. Making 

papu _ Ufl, 

■\x\ ~ P\x\ - P\x\ 



use of the symmetry, orthogonality, and tracelessness properties pf?" po ,r' 7 )"' , ■ 



P fJj p^ pcr = 0, and g pu p p ^ pc7 = 0, respectively, we can constrain the decomposition of S p,vpa 
to be of the form 
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gpupa 



pupa 



18) 



i=l 



where F and q are scalar functions of s, t, and u, defined in Eq. fllOD, and 



pupa 
S l 

pupa 
s 2 

pupa 
S 3 

pupa 



pupa 
S 5 

pupa 
S 6 

pupa 
s 7 

pupa 



pupa 
S 9 

pupa 
b 10 



tt v P a 

PbPbPbPb , 

P"aP U aP P aPl 
U v P a 
PbPbPbPa, 

P"aP U aP P bPl 

PtPbP'aPl 



(19) 



The dependence on the reference frame of the measurement enters through Eq. (Q) and 
resides in the factors s ipu p p ^ and s ilMUpcT p p ^ pa , which are independent of the partonic sub- 
process. These factors explicitly depend on the longitudinal-momentum fractions which 
the partons a and b receive from the hadrons from which they spring. Their calculation 
is straightforward, and we refrain from listing them here. The functions F and q in 
Eqs. (|^) and ( |TE| ) are independent of the hadron momenta and the reference frame; they 
only depend on the partonic subprocess. The set of these functions which are relevant for 
the direct hadroproduction of polarized 3 Si charmonium (H = J/ip,ip') in the NRQCD 
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factorization formalism may be found in Refs. PJTT|. This includes the 3 5j channel, 
which also contributes in the cases H = Xci, Xc2- In Appendix |A], we provide the missing 
sets of functions pertinent to the latter two cases, namely those associated with the 3 P X 
and 3 p2 channels, respectively. 



3 Polarization in Cascade 

We now explain how to calculate the polarization variable a, defined in Eq. ([!]), for prompt 
J /if) production in the NRQCD factorization formalism. Our approach not only applies to 
hadroproduction, but also to any other production mechanism, such as photoproduction, 
deep-inelastic scattering, e + e~ annihilation, two-photon scattering, etc. It is convenient 
to calculate a from the unpolarized cross section a and its longitudinal component o L . 
Prompt J /if) mesons originate from direct J /if} production, from single-cascade decays of 
directly produced Xcj and if)' mesons, and from double-cascade decays of directly pro- 
duced if>' mesons via Xcj intermediate states. The unpolarized cross section crP rom P t is 
simply obtained by adding the unpolarized cross sections of the various direct-production 
processes multiplied with the appropriate branching fractions. 

The calculation of the longitudinal cross section a p * ompt j s m ore involved. In the 
following, we explain how it can be expressed in terms of the reduced partonic cross 
sections <7|^|(n) for the production of cc pairs in Fock states n with definite absolute 
helicities |A|, scalar ME's, and branching fractions. We have 

^prompt J ft = ^direct J ft + ^ + ^ + ^XcJ (2Q) 

The cross sections for the direct production of polarized J/if> and if>' mesons is given by 



a 



direct tp 
L,T 



E*0»(O*(»)>. (21) 



where the summation is over n = 3 S[ 1 \ 3 S[ 8 \ 1 Sq 8 \ and 3 Pj 8 ' ) . The relevant reduced 
partonic cross sections <7o,i(n) may be extracted from Refs. IHPI- Equation (^T|) allows 
us to obtain the first term in Eq. (|20|) . 



In the NRQCD factorization formalism, the radiative decay Xcj — * J/if> + 7 is ap- 
proximately treated as an electric dipole transition. Due to heavy-quark spin symmetry, 
the helicity of the J /if) meson is then given by the third spin component ms of the Xcj 
meson. Thus, angular-momentum conservation constrains the transition probability of a 
Xcj meson with helicity mj to a longitudinally polarized J/if> meson, with ms = 0, to 
be proportional to Y^ln L =-i \ (^-, m L', 1 , 1 «/, mj)\ 2 . By applying this kind of analysis to the 
leading cc Fock states n = s Pj^ and 3 S^ relevant for the Xcj mesons, we find the second 
term in Eq. (|20| ) to be 



a x CJ 



\a ( 3 P (1) ) (O** ( 3 P (1) )) + \a ( 3 5i 8) ) (0- ( 3 si 8) 



7 



i 



i 



3o(8) 



(8) N 



x B(xci -> ■//> + 7) 
+ 



( 3 P 2 (1) J + (^ 2 



3pW 



o ft p (l) 



£)XcO ip 







' ^° W + ( 3 ^ 8) )1 C^)) - W + 7)- (22) 



.30 u v 1 V 60 
Here, we made use of the multiplicity relations 

(pxcj ( 3 P«)^ = (2 J + 1) (p Xc0 ( 3 P {1) 

ft^)) ^ = ( 2 J + 1) (flXo ft^ ^ ^ (23) 

which follow to leading order in v 2 from heavy-quark spin symmetry. To leading order 
in v 2 , the absolute value of the first derivative of the radial wave function of P-wave 
charmonium at the origin |Pp(0)|, which is one of the input parameters of the CSM, is 
related to (O Xc0 ( 3 P (1) )) by (C Xc0 ( 3 ^o (1) )) = 3JV C \R' P (0)\ /(2tt), where N c = 3. The 

j J agree 



combinatorial factors multiplying the terms in Eq. (|2~2|) which involve a\\\ ( 3 -P 
with those found in Ref. 



The hadronic decay ip' — ► J /if) proceeds predominantly through a double chromoelec- 
tric dipole transition. Since transitions involving a spin flip are suppressed and the recoil 
momentum of the J/ip meson is negligible, the J/ip meson has the same helicity as the 
ip' meson. Thus, the third term in Eq. (|20"D may be evaluated as 



a 



direct i/j' 



Bty' -> J/ip + X). 



(24) 



The fourth term in Eq. (f20|) may be derived in a similar fashion as the contribution 
to <j x i cJ from the 3 S , j 8 ^ channel, and it reads 

*V XcJ = 1^ *B{# - Xco + l)B( Xc0 - J/V + 7) 

+ (iof cct * + iof* P(^' - Xd + l)B( X ci - W + 7) 

+ ^T ircct W - Xc2 + l)B{ Xc 2 - J/4> + 7). (25) 



i 1' direct lb' 

'30^ 



4 Numerical Analysis 

We are now in a position to explore the phenomenological implications of the analytic 
results presented in Sees. ^| and ^|, and Appendix |A] for the inclusive hadroproduction 
of polarized J/ip mesons at the Tevatron, with CM energy yfs = 1.8 TeV. As in the 
CDF analysis ||, we consider the px distribution da/dpx, integrated over the rapidity 
interval \y\ < 0.6. We work at leading order in the NRQCD factorization formalism. We 



use the MRST98LO parton density functions (PDF's) of the proton [14j] and evaluate 
a s (fj,) from the one-loop formula with Aqcd = 174 MeV fllil . We set the renormalization 
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and factorization scales \i equal to the transverse mass mx = y4m;? +Pr> with m c = 
1.5 GeV. We adopt the nonperturbative ME's appropriate for our choice of PDF's from 
Table I of Ref. |J. As for the J/ip and ip 1 mesons, fits to the Tevatron data do not 
simultaneously constrain (O^ (^Sq 8 ^^ and ((D^ ( 3 -Po j)j but only the linear combination 
M r = (p* + r (p* ( 3 Pq S) )) /m% with some optimal value of r. For simplicity, 

we assume that (p* ( 1 ^ 8) )} = r (O^ ( 3 P (8) )) l m l = M r/2. We adopt the values of the 
relevant branching fractions from Ref. We essentially work within the fusion picture, 
where the cc bound state is formed in the primary hard-scattering process. In the high-py 
regime, only those partonic subprocesses survive where the cc pair is created from a single 
gluon which is close to its mass shell. These contributions, which contain large logarithms 
of the form ln(mr/2m c ), are conveniently implemented by using fragmentation functions 
(FF's). These logarithms are then resummed by evolving the FF's from their starting 
scale \i = 2m c up to the scale \i = mr- In the case of the unpolarized cross section 
^prompt j/i> we a pp rox j m ately include this resummation by multiplying the fusion cross 
section a P-S^J with the ratio of the corresponding fragmentation cross sections with 
final-state factorization scales /i = my and \x = 2m c . In the case of the longitudinal 
cross section ,jP rompt we include the LO fragmentation cross sections o"o(n) with n = 
3 S{ 8 \ 1 SP, and 3 Pj 8) @, which are of order a*. These contributions are particularly 
important because the corresponding fusion cross sections vanish to order a 3 . Further 
details concerning the implementation of the fragmentation contributions may be found 
in Ref. §. 

In Fig. |l|, we study the pt dependence of the polarization variable a for prompt J/ip 
mesons. For comparison, we also show a for the various production mechanisms which 
contribute to the prompt- J/ip signal, namely direct J /if) production, single-cascade decays 
of directly produced Xcj and ip' mesons, and double-cascade decays of directly produced ip' 
mesons via Xcj intermediate states. A quantitative understanding of the asymptotic large- 
Pt behaviour of a may be obtained by observing that, in this limit, inclusive charmonium 
production dominantly proceeds via the gluon fragmentation process g cc (^S^^j |T7] . 
Since large-p^ gluons are almost on shell and thus predominantly transversely polarized, 
a i ( 3 s[ 8 ^ should greatly exceed all other reduced partonic cross sections at large pt values 



fl8|| . It hence follows that, for directly produced J/ip mesons, a approaches unity in the 
large-pr limit || . By heavy-quark spin symmetry, the same is true for J / ip mesons from 
single-cascade decays of ip 1 mesons ||. The situation is more involved for J/ip mesons 
from the radiative feed-down of Xcj mesons and for those from double-cascade decays of ip' 
mesons. Inspecting Eqs. ( p2|) and ( Iffi) in the large-pr limit, we can derive the asymptotic 
expression 

_ 60a:o + 15xi + 21x 2 
01 ~ 60x + 75xi + 73x 2 ' 

where xj = (2J+l)B(xcj — > J/'P + l) in the first case and x j = B{ip' — > Xcj + l)B(XcJ - ► 
J/ip + 7) in the second one. Experimental data suggests that B(xco — > J/'P + 7) "C 
B(Xa ~* J/'P + 7) ~ 2-B(xc2 J/'P + 7) and that B{ip' — > Xcj + 7) is approximately 
independent of J WB. We thus obtain a x - J » 39/163 w 0.24 and a^'^ x - J w 51/223 » 
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0.23. (Without these approximations, we find a XcJ ~ 0.24 and a^'^ XcJ ~ 0.24.) From 
Fig. [I] we see that the corresponding predictions at px = 30 GeV are already very close 
to these asymptotic values. Finally, we observe that the predicted pt distribution of 
a from prompt J/ip mesons is relatively close to the prediction for direct J/ip mesons. 
As is evident from Fig. [1], this is a consequence of a strong cancellation between the 
Xcj J/i> + 7 and ip' — > J/ip + X channels, which could not have been anticipated 
without explicit calculation. 

5 Conclusions 

The very recent measurement of the polarization of prompt J/ip mesons at the Tevatron 
offers a unique opportunity to test the NRQCD factorization formalism [|l|, which 
postulates the existence of color-octet processes in nature. The theoretical prediction of 
prompt- J/ip polarization is complicated by the contributions due to the feed-down from 
heavier charmonium states, which come in addition to the direct- J/ip contribution and 
make up about 40% of the prompt- J/ip rate. Specifically, one needs to consider the single- 
cascade decays of directly produced Xci, Xc2, and ip' mesons, and the double- cascade decays 
of directly produced ip' mesons via Xcj intermediate states. In this paper, we provided, 
in analytic form, the missing ingredients which are necessary to make this prediction. 
These include the cross sections of the partonic subprocesses a + b — > cc(n) + d with 
n = and 3 i~2 , and the probabilities for the transitions from the cc Fock states of 

the Xci, Xc2, and ip' mesons to those of the J/ip meson in the presence of polarization. We 
presented quantitative predictions, in terms of the polarization variable a, for the various 
production mechanisms which contribute to the prompt- J/ip signal. In the large-pr limit, 
a approaches the asymptotic values a direct J ^ — oft = 1 and a XcJ ~ ^ XcJ ~ 0.24. 

A comprehensive discussion of a for prompt J / ip mesons together with a conservative 
analysis of the theoretical errors may be found in Ref. M. The theoretical prediction 
is consistent with the CDF data at intermediate values of pr, but it exceeds the 
data point of largest pt by about three standard deviations of the latter. For the time 
being, it is premature to conclude that this represents experimental evidence against 
the NRQCD factorization formalism. On the theoretical side, the prediction should be 
improved by the inclusion of next-to-leading-order corrections in a s and v 2 , which are 
not yet available. This could have an appreciable effect, since a s (2m c ) « v 2 ~ 0.3 is 
not actually so small against unity. On the experimental side, Run II of the Tevatron is 
expected to increase the presently available charmonium data sample by roughly a factor 
of 20. This will allow for the polarization to be measured with higher precision and out 
to larger values of px- Furthermore, it should then be feasible to discriminate between 
the various channels which constitute the prompt- J/ip signal, as is presently done in the 
measurement of the unpolarized cross section. Should the discrepancy persist despite 
concerted theoretical and experimental progress, then one possible explanation could be 
that the NRQCD factorization formalism fails for charmonium simply because m c is too 
small. Other effects, such as the one due to nonvanishing intrinsic transverse momenta of 
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the initial-state gluons might be important. It is, therefore, worthwhile to perform 
parallel studies for polarized bottomonium. 
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A Partonic Cross Sections for Color-Singlet Produc- 
tion of Polarized Xci an d Xc2 Mesons 

In this appendix, we present the ingredients which are necessary to describe the direct 
hadroproduction of polarized xa and Xc2 mesons via the color-singlet channels n = 3 P^ 
and 3 i 3 2 j respectively. The reduced differential cross sections dm\\/di of the partonic 
subprocesses a + b —>■ cc(n) + d, where the cc pair has absolute helicity |A|, are given by 

^ = AjM. (A.l) 

dt 167TS 2 ' 'l A l V ^ 

Here, the partons a, b, and d are treated as massless, while the charm quark c has mass 
m c . In the spin-one (spin-two) case, |A^|m is given by Eq. (§) (Eq. JT7|)), which contains 
the process-dependent tensor S^ u (S^ vpa ). The functions F and q which appear in the 
Lorentz decomposition of the latter according to Eq. (||) (Eq. flUf)) are listed below. As 
in Eq. (|2"2"D, the color-sing let ME (p*» ( 3 ^o (1) )) has been factored out. Summing over 
|A|, we recover the well-known results for the unpolarized case [Q. As a further check 
for our analytical analysis, we also verified the differential cross sections of the partonic 
subprocesses a + b — >• cc ( 3 Po ) relevant for the hadroproduction of Xco mesons, given 
in Ref. f|. 
qq — > cc f 3 -Pi ) + g: 

64(4vrcv s ) 3 



81m c (t + u) 4 
Ci = tu, 
c 2 = 2t, 
c 3 = 2u, 

c 4 = 2(t + w). (A.2) 
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cc 



cc^P^+q: 



cc( 3 P 2 (1) ) +q: 



cc^P^+g: 



F = 



128(47ra s ) 3 m c 



s 2 + i 2 + 2s(i + u) 
s 2 + ii 2 + 2s(i + u 



81s 2 (t + w) 4 ' 
ci = s(i + u) 2 , 

C 2 =4: 

c 3 = 4 l 

c 4 = 8 s 2 -tw + 2s(t + w) 
C5 = c 6 = 8s, 
c 7 = c 8 = 16s, 

Cg = — 16(£ + u), 

C10 = 16(s + i + u). 



F = 



8(47ra s ) ; 



27m c (s + w) 4 ' 
Ci = su, 
c 2 = 2s, 
c 3 = 0, 

C4 = 2(s — u). 

F _ lQ{A-na s fm c 
27t 2 (s + w) 4 ' 
ci = i(s + -u) 2 , 
c 2 = 4 [s 2 + t 2 + 2t(s + u) 

c 3 = 4s 2 m 2 , 
c 4 = 8s(s + -u), 
c 5 = 8f , 
c 6 = c 8 = 0, 
C7 = 16t, 

Cg = — 16(s + u), 

C10 = 16(s + i + u). 



6m c (s + t) 4 (t + u) A {u + s) 4 ' 
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(A.6) 



ci = -2s 2 i 2 u 2 (si + iu + us) (s 2 + i 2 + u 2 - 2si - 2tu - 2us) 

- s 2 Pu 2 (§ + i + u) \si(s + i) + iu(i + u) + us(u + s) 

- (si + iu + us) [s 3 t 3 (s - i) 2 + i 3 u 3 (i - u) 2 + u 3 s 3 (u - s) 2 ] , 
c 2 = 2(s + i) 2 [s 2 (i + u) [s(s - i)(s - u)(i + u) - iu(i - u) 2 

+ (i - u) [t 2 + u 2 — s(i + «)] [s 2 £ 2 + tV + w 2 s 2 + stw(s + 2t)] } , 
cs = c 2 \ { «u, 

c 4 = 4(t + u) 2 {s 3 (s 4 + t 4 + u 4 ) - siu [s 4 + s 2 (5t 2 - 7tu + 5m 2 ) + iu(i 2 - 3iu + u 2 

- 4(t + u) (i - u) 2 \s 4 (i + 2u) (2i + u)- s 2 iu(P + u 2 ) + £ 3 w 3 . 

99^cc( 3 P^)+g: 

F 8(47ra s ) 3 m c 

~ 3(siu) 2 (s + t) 4 (t + u) 4 (u + s) 4 ' 
ci = s£m(s + t) 2 (t + u) 2 (u + s) 2 (s 2 + 1 2 + u 2 )(si + iu + ws) 2 , 
c 2 = 2(s + i) 2 {Au 8 si(s + t) 2 + u 7 si(s + t) (7s + 26st + 7?) 

+w 5 s 2 t 2 (s + £)(32s 2 + 71st + 32?) + ifs 2 ? [lO(s 4 + i 4 ) + 64si(s 2 + i 2 ) + 105s 2 P 
+u 3 s 3 i 3 (s + t)(17s 2 + 32st + 17t 2 ) + w 2 s 4 t 4 (8s 2 + list + 8t 2 ) - us 5 i 5 (s + i) 
-s 6 t 6 } , 

c 3 = c 2 | t -^, 

c 4 = -2(s + £)(s + u) {s 6 (P - w 2 ) 2 (t + 2u) (2i + u) - s 5 t 2 w 2 (t + u)(3* + 5u)(5i + 3u) 
-s 4 tw [3(t 6 + u 6 ) + 28M(t 4 + u 4 ) + 83t 2 w 2 (t 2 + u 2 ) + 114t 3 M 3 " 
-s 3 ?m 2 [ll(£ 5 + u 5 ) + 53tu(£ 3 + u 3 ) + 92t 2 w 2 (t + u)" 

-s 2 i 3 u 3 (i + w) 2 (13t 2 + lOtu + 13m 2 ) - si 4 u 4 (t + w)(5t 2 + 6tu + 5m 2 ) - 2t 6 w 6 } , 
c 5 = 4stw 2 (s + i) 2 (2si + iu + us) 2 \2u(s + i + u) 2 + (s + i)(si - u 2 )} , 

c 7 = 4stw 2 (s + t) [t 5 (s + u)(5s 2 + 12su + 5m 2 ) + s 3 w 3 (3s 2 + Asu + 3u 2 ) 

+w 5 t(lls 2 + 13si + hi 2 ) + m 4 £(37s 3 + 68s 2 £ + 43s£ 2 + 8t 3 ) 

+w 3 st(43s 3 + 132s 2 t + 124st 2 + 45t 3 ) + w 2 s 2 t(19s 3 + 110s 2 t + 156s* 2 + 82t 3 ) 

+ms 3 P(31s 2 + 83s£ + 61?) + s 4 t 3 (lls + 16*)] , 

Cs = C7 1 

c 9 = 4(s + i)(s + u) {s 5 [2(t 6 + u 6 ) + 5iu(i 4 + u 4 ) + Ai 2 u 2 (i 2 + iu + u 2 ) 



+s 4 iu{i + u) [3(t 4 + m 4 ) - iu(i + u) 2 



2s 3 i 2 u 2 



i 4 + u 4 + 6iu(i 2 + u 2 ) + 9t 2 u 



1+2-2 
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f3 u 3 



-s 2 t 3 u 3 (t + u){t 2 + u 2 ) + Ast 4 u 4 {t + uf + 2t b u\t + u)}, 
c 10 = -4 {s 7 [2(t 6 + u 6 ) + 2tu{t 4 + u 4 ) - 2li 2 u 2 {t 2 + u 2 ) - 46t 

+s 6 (t + u) [2(£ 6 + m 6 ) + 4tw(t 4 + m 4 ) - 53t 2 w 2 (t 2 + w 2 ) - 134t 3 w 3 

2 ^ 428t 3 w 3 



+s 5 iu 



2(F 



u ' 



A3tu(t 4 



u 



256t 2 u 2 (t 2 



u 



-s 4 i 2 u 2 (t + u) [l5(t 4 + u 4 ) + imuit 1 + u 2 ) + 264t 2 u 2 " 

-s 3 i 3 u 3 (t + u) 2 (32i 2 + 93iu + 32u 2 ) - s 2 i 4 u 4 (i + n)(19t 2 + 34tu + 19m 2 



+2t 6 M 6 (s + t + w)}. 



(A.7) 
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Figure 1: Polarization variable a as a function of pt for direct J/ip mesons, for those from 
single-cascade decays of Xcj and ip' mesons, for those from double-cascade decays of ip' 
mesons via Xcj intermediate states, and for prompt ones. 
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